Galactic Streams of Cosmic-ray Electrons and Positrons 



Matthew D. Kistler, 1,2 Hasan Yiiksel, 3 and Alexander Friedland 3 

1 Lawrence Berkeley National Laboratory and Department of Physics, University of California, Berkeley, CA 94-720 

^Einstein Fellow 

3 Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87544 

(Dated: October 29, 2012) 

Isotropy is a key assumption in many models of cosmic-ray electrons and positrons. We find that 
simulation results imply a critical energy of ~ 10-1000 GeV above which e can spend their entire 
lives in streams threading magnetic fields, due to energy losses. This would restrict the number of 
e sources contributing at Earth, likely leading to smooth e spectra, as is observed. For positrons, 
this could be as few as one, with an enhanced flux that would ease energetics concerns of a pulsar 
origin of the positron excess, or even zero, bringing dark matter into play. We conclude that ideas 
about e propagation based on either isotropic diffusion or turbulent fields must be changed. 



PACS numbers: 95.85.Ry, 98.35.Eg, 98.70.-f, 98.70.Sa 

Introduction. — Energy losses due to synchrotron ra- 
diation and inverse-Compton scattering severely limit the 
distances over which cosmic-ray electrons and positrons 
can remain highly energetic [HH]. Thus, observations of 
these particles provide an opportunity to examine the re- 
cent high-energy history of within a few hundred parsecs 
of Earth. A number of experiments have contributed to 
the recent developments in this area [3HZ] and more are 
poised to break into the multi-TeV regime [8HT0l. The 
objects that may yield measurable fluxes in this range are 
few [11] and represent perhaps the most immediate hope 
of directly discovering an extrasolar cosmic-ray source. 

Making use of this data requires understanding of the 
propagation of cosmic rays through magnetic fields af- 
ter escaping their sources. Many standard treatments 
for electrons and positrons (e.g., Refs. [T2HT6"] ) assume 
isotropic diffusion. But why? In smooth fields, particles 
follow field lines, resulting in funneling in specific direc- 
tions. As we will see, considering small scale turbulent 
variation in the field |17j alone ends up being too weak to 
cause isotropization for physically sensible assumptions 
in scenarios involving e . 

We compare simulations of e ± propagation using ran- 
dom magnetic fields with such analytical solutions over a 
range of energies, distances, and times. We concentrate 
on particular realizations of fields, rather than the aver- 
age over many such fields as implicitly assumed in diffu- 
sion models, since we can observe particles directly from 
only one location in the Galaxy. As an illustration, we 
show in Fig.[l]particle trajectories arising from nine sim- 
ulated e ± sources distributed within one random field, 
as described in detail below. We explicitly see, similar 
to [18], that rather than being distributed in spherically 
symmetric diffusive distributions, the particles are ini- 
tially confined to filamentary structures, which we refer 
to as "streams" . 

Importantly, we take into account particle energy 
losses, which is crucial here, since cooling necessitates 
that electrons and positrons can only travel for a limited 
time with energies exceeding a given value, in stark con- 
trast to the case of protons. We argue that, because of the 



losses, the consequence of this approach is that electrons 
and positrons with measured energies > 10-1000 GeV can 
be expected to spend their entire lives within streams, 
never reaching the isotropic diffusive regime. 

This containment of cosmic-ray electrons and positrons 
would be far reaching, including a need to abandon 
isotropic e ± diffusion calculations in this regime. We 
discuss how this limits the number of sources that can 
contribute to the measured e spectra, which has the 
advantage of accounting for the absence of features ex- 
pected from the multitude of potential nearby sources. 
This would affect models for the positron excess, with 
implications for a dark matter origin, and we further ad- 
dress additional factors that should be examined. 
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FIG. 1: Distributions of electrons with initial energies of 
1 PeV after propagating 5000 yr in a 3 uG random magnetic 
field. Each stream (dots) corresponds to one of nine sources. 
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Fields and streams. — To study cosmic-ray electrons 
traversing a number of magnetic field realizations, we 
simulate the propagation of large numbers of particles to 
obtain particle densities over a range of times and dis- 
tances. The Galactic field consists of a regular compo- 
nent with reversals in direction between spiral arms and 
a turbulent component with coherence length ~ 100 pc 
[19j . The dominant field surrounding the Sun is tilted 
[20] . with interpretations of the IBEX Ribbon yielding 
a ~ 3 /iG field oriented ~ 55° off of the Galactic plane 
[2"Tl |2"2"] . Fluctuations consistent with Kolmogorov tur- 
bulence are inferred over a wide range of scales [23] . Ra- 
dio synchrotron data imply a local ~ 7.5 /iG random field 
[21] , and a few /iG turbulent field may explain small-scale 
PeV proton anisotropies [25] . 

For a singly-charged particle within a magnetic field of 
strength B, the equation of motion can be written as 
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where (3 is a unit vector since the velocity is given in 
terms of c. Both distance and time are expressed in pc, 
energy in PeV, and B in /iG, with a Larmor radius of 
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We first follow the formulation laid out in Ref. 
(see also [27H30] ) to obtain a 

between scales fco and k\ whose power follows 
£-(n+2) ( n = 5/3 f or Kolmogorov turbulence), address- 
ing other possibilities later. It is a practical impossibility 
to construct a magnetic field to the k values associated 
with the ri, of GeV-TeV particles, due to memory con- 
straints. We instead reach smaller scales by using nested 
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FIG. 2: Radial distributions of cosmic-ray e ± with initial en- 
ergies of 10 TeV-10 PeV after propagation through a 3/tG 
random magnetic field at given times (bins). Also shown are 
the expectations from isotropic diffusion (dashed) and Jiittner 
(solid) models, which are seen to poorly represent the data. 



boxes, B^WocEilo^Bfe^)- where V = k /h 
and boxes repeat periodically when needed (similar to 
the method used in Ref. [31]). This maintains the proper 
normalization of power over all scales. 

A sufficiently large ko value ensures variation in the 
simulation box even in the largest scales and rj can be as 
large as memory allows. We normalize the rms value 
of B 1 ^ to 3 /iG and choose a maximum length scale, 
Imax oc 1 /ko , of 200 pc and a minimum scale, l m in °c 1/fcjv , 
smaller than rx,, with a coherence length l c ~ 40 pc. Par- 
ticles are injected isotropically at the source and Eq. 
is solved via a fourth-order Runge-Kutta method with 
step sizes smaller by at least a factor of several than n,. 

Our interest is principally in high-energy electrons, 
for which energy losses arise from inverse-Compton (IC) 
scattering on ambient photons and synchrotron radia- 
tion. We approximate energy losses as continuous (IC 
is stochastic at very high energies |32j : however, gamma 
rays are not our focus), with —dE/dt — b(E) = b E 2 , so 
1/ 'E = 1/ 'Eg + bot, where E g is the energy at generation. 
Assuming a ~ 3 /iG magnetic field and the cosmic mi- 
crowave background (CMB), 6 ~5 x 10~ 17 s _1 GeV" 1 . 

In Fig. [l] we show the positions of cosmic-ray e ± with 
E g = 1 PeV after propagating 5000 yr from sources at with X(E, t) = J Q Z dt'D[E(t')} = f* 9 dE'D(E') /b(E') giving 



nine locations distributed within the same 3 /iG field con- 
figuration. We see that, rather than diffusing isotropi- 
cally, the fluxes tend to escape along specific directions. 
This has been seen in the context of sources of Galactic 
protons with similar parameters |18j . 

Using the source position located at the center of 
Fig. [TJ we consider instantaneous bursts for E g cover- 
ing 10 TeV-10 PeV (variable sources can be built up as 
combinations of successive bursts). We display in Fig. [2] 
snapshots of the radial distributions summed over all di- 
rections for these four energies at three given times. 

Comparison to analytical models. — Assuming 
spherical symmetry, the diffusion equation yields a par- 
ticle density nd(r,t,E) from a bursting source pQ as 
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FIG. 3: Electron streams arising after propagating for 5000 yr with initial energies of 1 PeV. Here, we vary the ratio of regular 
to random field magnitudes as (white), 1 (red), and 5 (blue), while fixing B TOS + -Brand = 3 /iG. We see that increasing B TBg 
orients propagation along the regular field direction (as indicated). 



^dif (E, t) — 2yJ X(E, t). Eg is mapped to the measured E 
after losses via dE g /dE~ (E g /E) 2 , using the above bo- 

Since the diffusive solution has been noted to fail at 
high energies [H] [33] , we also consider a model based on 
the Jiittner particle distribution, which explicitly con- 
strains v < c to eliminate superluminal behavior, giving 



to the diffusive regime for protons was found as 
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where £(r, t) — r/ct, 6 is the step function, K\ is the mod- 
ified Bessel function, and a,j(E, t) = c 2 t 2 /(2X(E, t)) [33]. 

In Fig. [2] we compare the simulated densities to both 
models, using D(E) = 5.4 x 10 26 (E/GeV) 1 ^ 3 cm 2 s"\ a 
value close to that seen in Ref. [IS]. We note that this 
is lower by about an order of magnitude than used in 
phenomenological diffusion models [33]. Neither provide 
a good representation of the radial distributions, nor can 
they convey the inherent asymmetries (as exhibited in 
Fig. [T I . The enhanced densities near the source position, 
resembling the expectation from diffusion in one dimen- 
sion, and concentrations in flux as compared to isotropic 
propagation lead us to describe these as streams. 

A short road to death. — We have seen, as in Fig. [2] 
that the appearance of streams remains pronounced as 
the particle energy is decreased. At sub-TeV energies, 
simulations become prohibitive, due to the decreasing e ± 
Larmor radius and the associated need for smaller time 
steps to accurately compute trajectories to determine the 
amount of time elapsed before propagation is isotropized. 
However, in Ref. [18) a scaling for the time of transition 
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where (3 ~ 2 and 7 = 0.25 — 0.5 for Kolmogorov (and 
similarly for other spectra). 

As already noted, energy losses are a vital consider- 
ation for electrons and positrons, since the energy loss 
timescale for high-energy e ± due to synchrotron radia- 
tion and inverse- Compton scattering is only 

<,~W^V^V(i^^lyr, (6) 
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where e 7 is the ambient photon energy density. It is 
apparent that lower energy than those considered in 
our simulations will be of interest. These do not ex- 
perience the Klein-Nishina suppression of higher-energy 
backgrounds, so we include estimates for the infrared and 
optical backgrounds near Earth [35 with the CMB, giv- 
ing e 7 ~ 1 eV cm -3 (neglecting variations with energy in 
the IC cross section). We see, for i max = 150-250 pc and 
Stat =4-7.5 nG, that tj = t d for E c 10 - 1000 GeV. 

The implication is that e with E > E c are expected 
to lose their energy prior to leaving streams and never 
reach the diffusive regime. This would have profound 
consequences for observing e ± with energies exceeding 
E c at Earth, since their densities are then governed by 
the structure of the local Galactic field. This is quite 
different than when dealing with protons, which experi- 
ence negligible energy loss and arrive from much greater 
distances after diffusing through many field domains. 
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What about B reg 1 — We have also examined the 
effect of including an overall regular field component. In 
Fig. [3j we show the outcomes from setting -B r eg/-B ran d 
i 1 . 1, and 5, while keeping B TCg + B ra nd — 3 fiG. We have 
again used the source positioned at the center of Fig. [T] 
with E g = l PeV. It is clear that, rather than diminish 
streams, including a regular field accentuates them. 

Results with a dominant random field, as for -Brand ~ 
7.5 fiG and B rcg ~ 3 /iG, are close to a pure random case. 
For -B reg ~ Brand, streams are aligned with neither the 
regular field nor the large-scale random components. In 
these scenarios, the direction of the nearby "regular" field 
likely lies off of the Galactic plane (as is observed). 

For larger B leg / B ran d, we see that propagation be- 
comes ever more one dimensional, with the stream con- 
tinuing to greater distances. We also find that the dif- 
ference in £)y compared to the phenomenological D van- 
ishes. Even when we decrease l max to 50 pc, for which 
Eqs. ([5} & ([6} imply a very large E c for a purely random 
field, we see the same basic behavior for B rog / B ran d = 
5 as in Fig. [3j Thus, in this regime the importance of 
streams would likely persist even at energies below the 
nominal E c associated with the random component. 

Discussion and conclusions. — For electrons and 
positrons at the energies of current and future interest, 
we have seen that a consideration of energy losses when 
propagation is treated as occurring in a turbulent mag- 
netic field alone can significantly change the expectations 
for observations, particularly for nearby sources. One 
consequence is a limitation to the number of contribut- 
ing sources at Earth. This holds when e propagation in 
the local Galactic magnetic field can be described using 
such a turbulent field coherent over a scale of ~ 100 pc. 
We discuss some of the implications that follow. 

Cosmic-ray electrons: As illustrated in Fig. [T] even for 
multiple local sources, at progressively higher energies 
one would be increasingly fortunate to be located in a 
stream. Crossed streams, where multiple sources would 
be inferred, are relatively rare. This is radically differ- 
ent than expected in isotropic models, where all nearby 
sources contribute to the flux and lead to a spectrum 
containing numerous features |16U36j . However, the mea- 
sured electron spectrum has a smooth variation [5J. 

These structures could be very important at the high- 
est energies for e , since energy losses are severe, and 
fluxes may be present even at TeV-PeV energies. Several 
nearby pulsars are known to accelerate e ± to these en- 
ergies [11] ; however, since streams do not follow straight 
lines, enhancements could appear in any direction, even 
when the source position is known. High-energy electron 
and positron anisotropics would be more complex than 
a simple dipole oriented with the source, as in diffusive 
prescriptions (see, e.g., [16l 137]). 

Positrons {and dark matter): Aside from the general 
rise associated with the positron excess 0] , positron data 
are also rather featureless and have been interpreted in 
terms of a single pulsar. Indeed, the smaller number of 



positron sources would imply fewer streams. One con- 
sequence could be an enhanced flux from an otherwise 
unremarkable source, easing concerns of exceeding the 
spin-down power of a lone pulsar. Limiting the sources 
that can contribute at Earth in this way could account 
for the smooth electron and positron spectra, since e ± 
fluxes from a source outside of a stream are diminished. 

The closer propagation is to isotropy, the more likely a 
particular source contributes a finite flux. We have seen 
in Fig. [T] that such "clouds" are occasionally realized. 
One can determine volume filling factors as a function 
of energy, distance, and time. However, one could argue 
that the positron excess suggests that at least one source 
does reach Earth, so we leave such probabilistic argu- 
ments to elsewhere to keep focus on the larger points. 

However, a scenario in which none of the local astro- 
physical positron sources actually reach Earth becomes 
a distinct possibility. At first glance, this would seem 
to require the diffuse injection of positrons, as from the 
smooth Galactic dark matter halo. However, positron 
models based on annihilations in dark matter substruc- 
ture (e.g., [38H4T] ) are affected, since each is a continu- 
ously emitting source. While it is even less likely for the 
largest "clumps" to contribute, the subhalo mass func- 
tion implies numerous less massive objects [I21II3]. Any 
of these could yield an apparent flux larger than expected 
from isotropic propagation (as assumed in, e.g., [4*4"l |4"5] ) . 

Our Galactic neighborhood: Streams would naturally 
limit the number of sources that contribute to the elec- 
tron and positron fluxes at Earth, without fine tuning 
the fluxes from myriad sources to achieve the smooth- 
ness present in the measured spectra. They could also 
ease energetics concerns by concentrating the flux from 
a source. However, we have seen this to be a significant 
departure from the usual picture, and since E c depends 
strongly on the actual parameters, which vary through- 
out the Galaxy, comparing to local data is essential. 

The complex nearby environment of Earth |46j can be 
uniquely probed by e due to their limited lifetimes. The 
above simulations contain the common assumption of a 
static field with no feedback effects. It is likely that 
propagation is more involved than this, and it is impor- 
tant to understand at what level. Depending on the in- 
terstellar medium, collective effects could increase scat- 
tering en route below some energy j47j [48] , as proposed 
to explain breaks in cosmic-ray protons [49 , or in the 
vicinity of sources [5"0"U5"!2] . We thus encourage their ex- 
amination in a consistent framework relevant near Earth 
that includes energy losses. In any case, we have found 
the investigation of e ± to be richer than often assumed 
and open to further development. 
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